Selection and Validation of Reference Genes in Sudan Grass (Sorghum sudanense (Piper) Stapf) under Various Abiotic Stresses by qRT-PCR.

Quantitative reverse transcription PCR (qRT-PCR) can screen applicable reference genes of species, and reference genes can be used to reduce experimental errors. Sudan grass (Sorghum sudanense (Piper) Stapf) is a high-yield, abiotic-tolerant annual high-quality forage with a wide range of uses. However, no studies have reported reference genes suitable for Sudan grass. Therefore, we found eight candidate reference genes, including UBQ10, HIS3, UBQ9, Isoform0012931, PP2A, ACP2, eIF4α, and Actin, under salt stress (NaCl), drought stress (DR), acid aluminum stress (AlCl3), and methyl jasmonate treatment (MeJA). By using geNorm, NormFinder, BestKeeper, and RefFinder, we ranked eight reference genes on the basis of their expression stabilities. The results indicated that the best reference gene was PP2A under all treatments. eIF4α can be used in CK, MeJA, NaCl, and DR. HIS3 can serve as the best reference gene in AlCl3. Two target genes (Isoform0007606 and Isoform0002387) belong to drought-stress-response genes, and they are highly expressed in Sudan grass according to transcriptome data. They were used to verify eight candidate reference genes under drought stress. The expression trends of the two most stable reference genes were similar, but the trend in expression for Actin showed a significant difference. The reference genes we screened provided valuable guidance for future research on Sudan grass.

Transcriptome Analysis Reveals the Potential Molecular Mechanisms of Tiller Bud Development in Orchardgrass.

Tillering is a special type of branching and one of the important contributors to the yield of cereal crops. Strigolactone and sucrose play a vital role in controlling tiller formation, but their mechanism has not been elucidated completely in most crops. Orchardgrass (Dactylis glomerata L.) is an important perennial forage with prominent tillering ability among crops. To date, the mechanism of tillering in orchardgrass is still largely unknown. Therefore, we performed a transcriptome and miRNA analysis to reveal the potential RNA mechanism of tiller formation under strigolactone and sucrose treatment in orchardgrass. Our results found that D3, COL5, NCED1, HXK7, miRNA4393-z, and miRNA531-z could be key factors to control tiller bud development in orchardgrass. In addition, strigolactones might affect the ABA biosynthesis pathway to regulate the tiller bud development of orchardgrass, which may be related to the expression changes in miRNA4393-z, NCED1, and D10. miRNA531-z could be involved in the interaction of strigolactones and sucrose in regulating tillering. These results will be further used to clarify the potential mechanism of tillering for breeding new high-tillering and high-production orchardgrass varieties and beneficial to improving the production and reproduction of crops.

Overexpression of Lolium multiflorum LmMYB1 Enhances Drought Tolerance in Transgenic Arabidopsis.

Lolium multiflorum is one of the world-famous forage grasses with rich biomass, fast growth rate and good nutritional quality. However, its growth and forage yield are often affected by drought, which is a major natural disaster all over the world. MYB transcription factors have some specific roles in response to drought stress, such as regulation of stomatal development and density, control of cell wall and root development. However, the biological function of MYB in L. multiflorum remains unclear. Previously, we elucidated the role of LmMYB1 in enhancing osmotic stress resistance in Saccharomyces cerevisiae. Here, this study elucidates the biological function of LmMYB1 in enhancing plant drought tolerance through an ABA-dependent pathway involving the regulation of cell wall development and stomatal density. After drought stress and ABA stress, the expression of LmMYB1 in L. multiflorum was significantly increased. Overexpression of LmMYB1 increased the survival rate of Arabidopsis thaliana under drought stress. Under drought conditions, expression levels of drought-responsive genes such as AtRD22, AtRAB and AtAREB were up-regulated in OE compared with those in WT. Further observation showed that the stomatal density of OE was reduced, which was associated with the up-regulated expression of cell wall-related pathway genes in the RNA-Seq results. In conclusion, this study confirmed the biological function of LmMYB1 in improving drought tolerance by mediating cell wall development through the ABA-dependent pathway and thereby affecting stomatal density.

Molecular Mechanism Underlying the Sorghum sudanense (Piper) Stapf. Response to Osmotic Stress Determined via Single-Molecule Real-Time Sequencing and Next-Generation Sequencing.

Drought, as a widespread environmental factor in nature, has become one of the most critical factors restricting the yield of forage grass. Sudangrass (Sorghum sudanense (Piper) Stapf.), as a tall and large grass, has a large biomass and is widely used as forage and biofuel. However, its growth and development are limited by drought stress. To obtain novel insight into the molecular mechanisms underlying the drought response and excavate drought tolerance genes in sudangrass, the first full-length transcriptome database of sudangrass under drought stress at different time points was constructed by combining single-molecule real-time sequencing (SMRT) and next-generation transcriptome sequencing (NGS). A total of 32.3 Gb of raw data was obtained, including 20,199 full-length transcripts with an average length of 1628 bp after assembly and correction. In total, 11,921 and 8559 up- and down-regulated differentially expressed genes were identified between the control group and plants subjected to drought stress. Additionally, 951 transcription factors belonging to 50 families and 358 alternative splicing events were found. A KEGG analysis of 158 core genes exhibiting continuous changes over time revealed that 'galactose metabolism' is a hub pathway and raffinose synthase 2 and ß-fructofuranosidase are key genes in the response to drought stress. This study revealed the molecular mechanism underlying drought tolerance in sudangrass. Furthermore, the genes identified in this study provide valuable resources for further research into the response to drought stress.

Integrated Analysis of Single-Molecule Real-Time Sequencing and Next-Generation Sequencing Eveals Insights into Drought Tolerance Mechanism of Lolium multiflorum.

Lolium multiflorum is widely planted in temperate and subtropical regions globally, and it has high economic value owing to its use as forage grass for a wide variety of livestock and poultry. However, drought seriously restricts its yield and quality. At present, owing to the lack of available genomic resources, many types of basic research cannot be conducted, which severely limits the in-depth functional analysis of genes in L. multiflorum. Therefore, we used single-molecule real-time (SMRT) and next-generation sequencing (NGS) to sequence the complex transcriptome of L. multiflorum under drought. We identified 41,141 DEGs in leaves, 35,559 DEGs in roots, respectively. Moreover, we identified 1243 alternative splicing events under drought. LmPIP5K9 produced two different transcripts with opposite expression patterns, possibly through the phospholipid signaling pathway or the negatively regulated sugar-mediated root growth response to drought stress, respectively. Additionally, 13,079 transcription factors in 90 families were obtained. An in-depth analysis of R2R3-MYB gene family members was performed to preliminarily demonstrate their functions by utilizing subcellular localization and overexpression in yeast. Our data make a significant contribution to the genetics of L. multiflorum, offering a current understanding of plant adaptation to drought stress.

Genome-wide identification of C2H2-type zinc finger gene family members and their expression during abiotic stress responses in orchardgrass (Dactylis glomerata).

The C2H2-type zinc finger protein (ZFP) family is one of the largest transcription factor families in the plant kingdom and its members are involved in plant growth, development, and stress responses. As an economically valuable perennial graminaceous forage crop, orchardgrass (Dactylis glomerata) is an important feedstuff resource owing to its high yield and quality. In this study, 125 C2H2-type ZFPs in orchardgrass (Dg-ZFPs) were identified and further classified by phylogenetic analysis. The members with similar gene structures were generally clustered into the same groups, with proteins containing the conserved QALGGH motif being concentrated in groups VIII and IX. Gene ontology and miRNA target analyses indicated that Dg-ZFPs likely perform diverse biological functions through their gene interactions. The RNA-seq data revealed differentially expressed genes across tissues and development phases, suggesting that some Dg-ZFPs might participate in growth and development regulation. Abiotic stress responses of Dg-ZFP genes were verified by qPCR and Saccharomyces cerevisiae transformation, revealing that Dg-ZFP125 could enhance the tolerance of yeasts to osmotic and salt stresses. Our study performed a novel systematic analysis of Dg-ZFPs in orchardgrass, providing a reference for this gene family in other grasses and revealing new insights for enhancing gene utilization.

MeJA-responsive bHLH transcription factor LjbHLH7 regulates cyanogenic glucoside biosynthesis in Lotus japonicus.

Cyanogenic glucosides (CNglcs) play an important role in plant defense response; however, the mechanism of regulation of CNglc synthesis by the external environment and endogenous hormones is largely unclear. In this study, we found that jasmonates (JAs) promoted the synthesis of CNglcs by activating the expression of CNglc biosynthesis genes in Lotus japonicus. Several differentially expressed basic helix-loop-helix (bHLH) family genes related to the synthesis of CNglcs were identified by RNA-seq. LjbHLH7 can directly activate the expression of CYP79D3 gene, the first step of CNglc synthesis, by binding to the G-box sequence of its promoter. Transgenic plants overexpressing LjbHLH7 exhibited higher relative CNglc content and enhanced insect resistance compared with the wild type. Furthermore, the transcriptional activity of LjbHLH7 was suppressed by the interaction with the L. japonicus JASMONATE-ZIM DOMAIN protein LjJAZ4. Based on these results, we propose that LjbHLH7 acts as an activator and LjJAZ4 acts as a repressor of JA-induced regulation of CNglc biosynthesis in L. japonicus.

Identification of novel drought-responsive miRNA regulatory network of drought stress response in common vetch (Vicia sativa).

Drought is among the most important natural disasters with severe effects on animals and plants. MicroRNAs are a class of noncoding RNAs that play a crucial role in plant growth, development, and response to stress factors, including drought. However, the microRNAs in drought responses in common vetch (Vicia sativa), an annual herbaceous leguminous plant commonly used for forage by including it in mixed seeding during winter and spring, have not been characterized. To explore the microRNAs' response to drought in common vetch, we sequenced 10 small RNA (sRNA) libraries by the next-generation sequencing technology. We obtained 379 known miRNAs belonging to 38 families and 47 novel miRNAs. The two groups had varying numbers of differentially expressed miRNAs: 85 in the comparison group D5 vs C5 and 38 in the comparison group D3 vs C3. Combined analysis of mRNA and miRNA in the same samples under drought treatment identified 318 different target genes of 123 miRNAs. Functional annotation of the target genes revealed that the miRNAs regulate drought-responsive genes, such as leucine-rich repeat receptor-like kinase-encoding genes (LRR-RLKs), ABC transporter G family member 1 (ABCG1), and MAG2-interacting protein 2 (MIP2). The genes were involved in various pathways, including cell wall biosynthesis, reactive oxygen removal, and protein transport. The findings in this study provide new insights into the miRNA-mediated regulatory networks of drought stress response in common vetch.

Genome-wide identification, phylogenetic analysis, and expression analysis of the SPL gene family in orchardgrass (Dactylis glomerata L.).

SPL (SQUAMOSA promoter binding protein-like) is a plant-specific transcription factor family that contains the conserved SBP domain, which plays a vital role in the vegetative-to-reproductive phase transition, flowering development and regulation, tillering/branching, and stress responses. Although the SPL family has been identified and characterized in various plant species, limited information about it has been obtained in orchardgrass, which is a critical forage crop worldwide. In this study, 17 putative DgSPL genes were identified among seven chromosomes, and seven groups that share similar gene structures and conserved motifs were determined by phylogenetic analysis. Of these, eight genes have potential target sites for miR156. cis-Element and gene ontology annotation analysis indicated DgSPLs may be involved in regulating development and abiotic stress responses. The expression patterns of eight DgSPL genes at five developmental stages, in five tissues, and under three stress conditions were determined by RNA-seq and qRT-PCR. These assays indicated DgSPLs are involved in vegetative-to-reproductive phase transition, floral development, and stress responses. The transient expression analysis in tobacco and heterologous expression assays in yeast indicated that miR156-targeted DG1G01828.1 and DG0G01071.1 are nucleus-localized proteins, that may respond to drought, salt, and heat stress. Our study represents the first systematic analysis of the SPL family in orchardgrass. This research provides a comprehensive assessment of the DgSPL family, which lays the foundation for further examination of the role of miR156/DgSPL in regulating development and stress responses in forages grasses.

Comparative transcriptome analyses reveal different mechanism of high- and low-tillering genotypes controlling tiller growth in orchardgrass (Dactylis glomerata L.).

BACKGROUND: Tillering is an important agronomic trait underlying the yields and reproduction of orchardgrass (Dactylis glomerata), an important perennial forage grass. Although some genes affecting tiller initiation have been identified, the tillering regulatory network is still largely unknown, especially in perennial forage grasses. Thus, unraveling the regulatory mechanisms of tillering in orchardgrass could be helpful in developing selective strategies for high-yield perennial grasses. In this study, we generated high-throughput RNA-sequencing data from multiple tissues of tillering stage plants to identify differentially expressed genes (DEGs) between high- and low-tillering orchardgrass genotypes. Gene Ontology and pathway enrichment analyses connecting the DEGs to tillering number diversity were conducted. RESULTS: In the present study, approximately 26,282 DEGs were identified between two orchardgrass genotypes, AKZ-NRGR667 (a high-tillering genotype) and D20170203 (a low-tillering genotype), which significantly differed in tiller number. Pathway enrichment analysis indicated that DEGs related to the biosynthesis of three classes of phytohormones, i.e., strigolactones (SLs), abscisic acid (ABA), and gibberellic acid (GA), as well as nitrogen metabolism dominated such differences between the high- and low-tillering genotypes. We also confirmed that under phosphorus deficiency, the expression level of the major SL biosynthesis genes encoding DWARF27 (D27), 9-cis-beta-carotene 9',10'-cleaving dioxygenase (CCD7), carlactone synthase (CCD8), and more axillary branching1 (MAX1) proteins in the high-tillering orchardgrass genotype increased more slowly relative to the low-tillering genotype. CONCLUSIONS: Here, we used transcriptomic data to study the tillering mechanism of perennial forage grasses. We demonstrated that differential expression patterns of genes involved in SL, ABA, and GA biosynthesis may differentiate high- and low-tillering orchardgrass genotypes at the tillering stage. Furthermore, the core SL biosynthesis-associated genes in high-tillering orchardgrass were more insensitive than the low-tillering genotype to phosphorus deficiency which can lead to increases in SL biosynthesis, raising the possibility that there may be distinct SL biosynthesis way in tillering regulation in orchardgrass. Our research has revealed some candidate genes involved in the regulation of tillering in perennial grasses that is available for establishment of new breeding resources for high-yield perennial grasses and will serve as a new resource for future studies into molecular mechanism of tillering regulation in orchardgrass.

Identification of candidate reference genes for quantitative RT-PCR in Miscanthus sinensis subjected to various abiotic stresses.

Quantitative real-time PCR (qRT-PCR) has been widely used for studying gene expression at the transcript level. Its accuracy usually relies on the reference genes that are utilized for data normalization. Miscanthus sinensis, a perennial C4 grass with high biomass and strong resistance to adversities, is often utilized as a high value energy crop. However, no reliable reference genes have been investigated for normalizing gene expression for this species. In this study, 12 candidate reference genes were selected to identify their stability under five different abiotic stress treatments (drought, salt, cadmium, chromium and arsenic) by using geNorm, NormFinder, BestKeeper and RefFinder softwares. The results showed that 18S rRNA and Unigene33312 were the best reference genes under drought treatments. Unigene33312 and Unigene33024 were found to be the most stably expressed genes under salt stress and Cd stress. Moreover, Unigene33024 and PP2A were the most suitable reference genes under Cr stress and Unigene33024 and Sb09g019750 were deemed more suitable reference genes under As stress. In total, considering all the samples, Unigene33024 and PP2A were the most stable genes while ACTIN and Unigene26576 were the least stable reference genes for internal control. The expression patterns of two target genes (Cu/Zn SOD and CAT) were used to further verify those selected reference genes under different conditions. The results showed that the most and the least stable reference genes had clearly different expression patterns. This work comprehensively estimated the stability of reference genes in M. sinensis which may give insight to the reference genes selection in other tissues as well as other related varieties. These suggested reference genes would assist in further putative gene expression validation in M. sinensis.

Genome-wide identification, structural analysis and expression profiles of GRAS gene family in orchardgrass.

The GRAS gene family is a family of transcription factors that regulates plant growth and development. Despite being well-studied in many plant species, little is known about this gene family in orchardgrass (Dactylis glomerata L.), one of the top four economically important perennial forage grasses cultivated worldwide. We identified 46 GRAS genes in orchardgrass and analyzed their characteristics by phylogenetic, gene structural, motifs and expression patterns analysis. The phylogenetic analysis of eight species revealed that DgGRAS family had the evolutional conservation and closer homology relationship with the GRAS family of rice, barley and Brachypodium distachyon. Moreover, 46 DgGRAS proteins were divided into eight subfamilies based on the tree topology and rice or Arabidopsis classification, and LISCL subfamily was the largest one. Besides, we found that the motif 15 may be unique to the orchardgrass LISCL subfamily, and the motif 6 and motif 17 had indispensable functions in the orchardgrass LISCL subfamily. We further analyzed the expression profiles of DgGRAS genes at mature and seeding stage. And we found that DgGRAS17 played an important role in the growth and development no matter what stage it was at. DgGRAS5, DgGRAS28, DgGRAS31, DgGRAS42 and DgGRAS44 got involved in processes of the growth and development at seeding stage instead of mature stage. These results indicated that the major expression patterns and detailed functions of the DgGRAS genes varied with developmental stages. Taken together, this is the first systematic analysis of the GRAS gene family in the orchardgrass genome and the results provide insights into the potential functions of GRAS genes.

De Novo Assembly and Discovery of Genes That Involved in Drought Tolerance in the Common Vetch.

The common vetch (Vicia sativa) is often used as feed for livestock because of its high nutritional value. However, drought stress reduces forage production through plant damage. Here, we studied the transcriptional profiles of common vetch exposed to drought in order to understand the molecular mechanisms of drought tolerance in this species. The genome of the common vetch has not been sequenced, therefore we used Illumina sequencing to generate de novo transcriptomes. Nearly 500 million clean reads were used to generate 174,636 transcripts, including 122,299 unigenes. In addition, 5313 transcription factors were identified and these transcription factors were classified into 79 different gene families. We also identified 11,181 SSR loci from di- to hexa-nucleotides whose repeat number was greater than five. On the basis of differentially expressed genes, Gene Ontology analysis identified many drought-relevant categories, including "oxidation-reduction process", "lipid metabolic process" and "oxidoreductase activity". In addition to these, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis identified pathways, such as "Plant hormone signal transduction", "Glycolysis/Gluconeogenesis" and "Phenylpropanoid biosynthesis", as differentially expressed in the plants exposed to drought. The expression results in this study will be useful for further extending our knowledge on the drought tolerance of common vetch.

Reference Gene Selection for Quantitative Real-Time Reverse-Transcriptase PCR in Annual Ryegrass (Lolium multiflorum) Subjected to Various Abiotic Stresses.

To select the most stable reference genes in annual ryegrass (Lolium multiflorum), we studied annual ryegrass leaf tissues exposed to various abiotic stresses by qRT-PCR and selected 11 candidate reference genes, i.e., 18S rRNA, E2, GAPDH, eIF4A, HIS3, SAMDC, TBP-1, Unigene71, Unigene77, Unigene755, and Unigene14912. We then used GeNorm, NormFinder, and BestKeeper to analyze the expression stability of these 11 genes, and used RefFinder to comprehensively rank genes according to stability. Under different stress conditions, the most suitable reference genes for studies of leaf tissues of annual ryegrass were different. The expression of the eIF4A gene was the most stable under drought stress. Under saline-alkali stress, Unigene14912 has the highest expression stability. Under acidic aluminum stress, SAMDC expression stability was highest. Under heavy metal stress, Unigene71 expression had the highest stability. According to the software analyses, Unigene14912, HIS3, and eIF4A were the most suitable for analyses of abiotic stress in tissues of annual ryegrass. GAPDH was the least suitable reference gene. In conclusion, selecting appropriate reference genes under abiotic stress not only improves the accuracy of annual ryegrass gene expression analyses, but also provides a theoretical reference for the development of reference genes in plants of the genus Lolium.